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Abstract 

The Inverse problem for an electromagnetic field produced by a dipole 
is solved. It is assumed that the field of an arbitrary changing dipole is 
known. Obtained formulae allow calculation of the position and dynamics 
of the dipole which produces the measured field. The derived results can 
be used in investigations on radiative process in solids caused by changing 
of the charge distribution. For example, generation of the electromagnetic 
field caused by oscillations of atoms or electron gas at the trace of a particle 
channeling in a crystal, or fields arising at solids cracking or dislocation 
formation - in any case when one is interested in the details of the dipole 
field source. 
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There are many applications where the problem of calculation of the source of 
the known electromagnetic field is arisen. For example, when the first undulators 
were developed, a question occurred on how to construct the undulator field in 
order to produce the radiation of given properties. Some ways of solution of this 
inverse problem was suggested in Refs m- An alternative approach concerning 
radiation in a ’’short magnet” was suggested by authors of Refs nig. The 
general solution of the inverse radiation problem is presented in Ref. [B]. The 
listed above papers present the methods of calculation of the particle motion 
if the produced electromagnetic field is known. The same problem is actual in 
astrophysics in research of distant sources of radiation. A general solution of 
the inverse problem for an arbitrary moving charged particle is given in Ref [7] . 

Another elementary and widespread source of electromagnetic field is a vari¬ 
able electric or magnetic dipole. Relativistic particle physics presents a great 
variety of the dipole radiation phenomena. One of the interesting example for 
inverse problem application is the field generated by atoms and electron gas 
exited by channeling particles in a crystal [5] . Of spatial interest are oscillations 
of atoms at a trace of a channeling particle. Solution of inverse problem for a 
field produced by an arbitrary changing dipole could be a good instrument of 
investigations in this area. Though the solution of the inverse problem for a 
dipole field is of great interest in very different areas of physics - for example, 
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in forecast of earthquakes or investigation of development of cracks in crystals 
(see a review in Ref. m) - the complete solution of the inverse problem for the 
dipole field has not been published until now. The solution in a particular case 
of the static dipole field was obtained in Refs ( [Min]). 

Electric E{t) and magnetic H{t) fields of a dipole momentum p{t') at a 
distance r is represented by the formulae: 
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where n = r/r, r is the radius-vector from the dipole moment to the observation 
point, c is the speed of light. Dots denote the derivatives with respect to time, 
the retarded time t' is connected with time t by relation 


t — t -\- 


r 

c 


(3) 


Eq. ([5]) shows that vector H is changing in a plane orthogonal to vector n, 
which allows, basically, to define the vector n. Eor this purpose we find the 
vector production {H x H). Taking derivative from Eq. ([5]) with respect to t 
and having in mind that dt = dt', we obtain: 

(H xH) = ^n{n ■ {P' x P")), (4) 

r‘ 


where P = p' + p, stroke denotes the derivative with respect to dimensionless 
time T = ctir. The left-hand side of the last equation gives the direction of 
angular velocity of vector E[. We see that both vectors El and P' rotate around 
the vector n in the same sense. The vector product {H xH) can be either parallel 
or antiparallel to the vector n depending on dipole dynamics. We assume further 
that {H X H) 7 ^ 0. Then we have for the unit vector n: 


ft = ± 


{H X H) 

\{Hxk)\ 


(5) 


The sign in this equation can be defined from the condition that the average 
energy flow defined by Poynting vector 

S = ±{ExH) (6) 

is directed from the dipole to the observation point. Indeed, let us find projection 
of S onto vector ft by use of Eqs CD and ([5]): 

{ft ■ S) = ^ g X p")^ + {ft X p') ■ {ft X p") -I- (n X P) • (n X P')| . 

The last equation can be written in the form 


{ft-S) 
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We assume that the dipole moment p varies within a limited interval. Then 
after avereging over time t from — oo to oo the value of second term terns to 
zero and the average rate of the first term represents the well known value of 
the dipole radiation: 


Thus, the true direction of vector n in Eq. is defined by condition {n ■ S) > 0. 
The second term in Eq. © can be both positive and negative, and describes 
the oscillatory part of the energy transition in the near-held zone. Generally 
speaking, direction of this part of energy flow is not parallel with the vector n. 

Now, knowing the vector n, we can integrate Eqs o and in order to 
calculate the dipole moment p{t'). After Hrst integration of Eq. ([2]) we obtain: 

#1 = 1(P X n), (8) 


where 


Hi 


/ 


H(t)dT 


is a function known to an additive constant. We adopt a Cartesian system of 
coordinates with the axis OX parallel to the vector n. Then Eq. (|5]) consist of 
two scalar equations: 


Hly — ^3 [Pz + Pz)i (9) 

Hiz = -^{p'y+Py)- (10) 

The general solution of this system of differential equations is as following: 

pj_ = Jin x Hi)e^dT^ , (11) 

where p± is the projection of vector p onto the plane YZ and poj_ is an arbitrary 
vector in the YZ plane. 

In order to hnd the X component of vector p we take a scalar production of 
Eq. (P) and vector ft. This gives: 

(”-^) = ^(P'x+P^)- 

Hence: 



where pox is an arbitrary constant. Combining Eqs (ED and ED we obtain the 
general solution of Eqs © and @: 

pit') = e~^ |po + J [in X Hi) + in(n • E)] e’^dr^ . (13) 
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Here po = po± + npo^ is an arbitrary vector. Up to now we have not calculated 
the distance r between the observation point and dipole. Hence we do not know 
the relation between times t in the right-hand side of Eq. m and t' in its 
left-hand side. 

Substituting the received solution (USD in initial Eqs UD and m we see that 
Eq. (I2D is satisfied identically but Eq. UD gives: 

E = n{n ■ E) — {n X H) -^e“'^po-L — J{n x Hi)e^dr. (14) 

The last relation shows that the vectors E and H are not independent - they 
are bound by Eq. USD- The relation between electric and magnetic fields of 
dipole can be found easily from Eqs UD and (ED by inspection. To do this we 
take a vector product of the vector n with Eq. ([TD and extract Eq. Q. As a 
result we get: 

H — {n X E) = x p). (15) 

Substituting (nxp) from this equation into Eq. we find the desired coupling 
equation: 


'^H+-H + H=—{nxE)E-{nxE). (16) 

c c c 

It is easy to show that Eq. (1141) can be transformed to Eq. ((T51) by multiplying 
into and taking the second derivative with respect to t. Thus, the solution 
USD satisfies the initial Eqs UD and ED- 

The distance between the dipole and observation point can be calculated by 
use of Eq. (ITSl) . It is a quadratic equation with respect to r: 


Ap^ + Bp + C = 0, (17) 

where vectors A = H — {n x E), B = H — {n x E), C = H are known, and 
p = r jc. The vectors A, B and C lie in a plane which is orthogonal to vector n. 
Hence, Eq. U3 is equal to the system of two scalar equations. Compatibility 
of the system requires certain relation between the coefficients of the equations 
which can be written in the form (A x C)^ = ((A x B) ■ {B x C)). Then Eq. 
UZD has the solution 

, An-jCxA)) ^n-jBxC)) 

(n-(AxH)) (n-(CxA)) 

In order to verify the obtained solution we represent the vectors A, B and C in 
terms of dipole moment. This can be made by use of Eq. UID: A = {px n)/r^, 
B = {p X n)/r^. The vector products take the form: 

_*_,77 ... T) ... 77 ... 

(AxH) = —{n-{pxp)), {CxA) = —{n-{pxp)), {BxC) = -^(n-{pxp)). 

The last relations show that the obtained solution for r is a positive constant, 
though the nominators and denominators in Eq. (jlSp are time-dependent. But 
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the expressions for r are not valid if (n ■ (p x p)) = 0, i.e. if vector p is varying 
in a plane going through vector n. In this case vectors P' P" are varying in 
the same plane, hence, the vector product in Eq. are equal to zero. As it 
was mentioned above this case is not under consideration in the current paper. 
Thus, the solution of the stated problem is given by Eqs (0, (HU) and m if 

{H X H) 7 ^ 0. And even if {H x H) = 0, Eq. (ITSl) for p(t') still remains valid. 

As an example we consider an electromagnetic field of oscillating dipole, say 
the field of an ion oscillating at specific site in the lattice of a crystal. Let the 
electric dipole have the time dependence 

Px = Pox coswt, py = poy sinut, Pz = Poz, (19) 

where pox, Poy, Poz, w are some given constants. We have calculated the electric 
and magnetic fields E{t) and H{t) at a distance r from the dipole and have saved 
the results as a numerical functions of t. Then we have used the saved data for 
calculation of p{t) and r{t) according to Eqs 0, (ITlTll and m- The result 
coincides with the initial functions with accuracy up to 10“®. The error occurs 
due to the used numerical method of calculation. The distance r due to the 
errors becomes a random function of time as shown in Fig. 1 for initial distance 
ruj/c = 100. 



Figure 1: The reconstructed distance r in terms of c/uj as a function of Time = 
cut 
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